Superfluid drag of two-species Bose-Einstein condensates in optical lattices 
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We study two-species Bose-Einstein condensates in quasi two-dimensional optical lattices of vary- 
ing geometry and potential depth. Based on the numerically exact Bloch and Wannier functions 
obtained using the plane-wave expansion method, we quantify the drag (entrainment coupling) 
between the condensate components. This drag originates from the (short range) interspecies in- 
teraction and increases with the kinetic energy. As a result of the interplay between interaction 
and kinetic energy effects, the superfiuid-drag coefficient shows a non-monotonic dependence on the 
lattice depth. To make contact with future experiments, we quantitatively investigate the drag for 
mass ratios corresponding to relevant atomic species. 

PACS numbers: 67.85.Fg, 67.85.Hj, 67.85.De 



I. INTRODUCTION 

Owing to their many favorable features such as the 
absence of defects and dynamic control of relevant pa- 
rameters, optical lattices constitute a versatile tool for 
controlling the properties of atomic quantum fluids [![. 
They have proven to be an invaluable platform for study- 
ing quantum phase transitions such as from a superfluid 
to a Mott-insulator (SF-MI) @, S] and simulating many- 
body systems jj| . The high densities and low tempera- 
tures required for investigating these quantum effects are 
provided by Bose-Einstein condensates (BECs) 0-0] ■ 

While multi-com pon ent BECs - and, especially, two- 
component ones [§-llfjj| - have attracted a great deal of 
interest since the early days of BEC jlj]: 13 L the addition 
of the optical-lattice environment [3|, [l4], [l5| allows one 
to study interesting transport properties of cold atoms 
in the superfluid regime. One of them, which has so 
far not been given due attention, is the superfluid drag 
in two-component BECs. Such a non-dissipative drag ef- 
fect was first investigated by Andreev and Bashkin [16[ in 
the context of 3 He- 4 He mixtures. A microscopic theory 
of the drag between two weakly interacting Bose gases 
was developed in the continuum limit [l7j and general- 
ized to the non Galilean invariant case as realized in an 
optical lattice [l8| ■ For a system of strongly interacting 
bosons, the superfluid drag was investigated by means of 
Monte Carlo simulations 19] . Drag effects have also been 
studied in different systems. For instance, in electronic 
mesoscopic systems Coulomb dra g ha s been studied both 
theoretically and experimentally [20| . 

In this work, we study the drag (entrainment coupling 
[2lj ) between two weakly interacting Bose gases in quasi 
two-dimensional optical lattices. We derive a formula for 
the supcrfluid-drag coefficient valid for an arbitrary lat- 
tice and evaluate it for different lattice geometries. In 
this manner we extend the results presented in Ref. [l8l 
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where this has been done for the special case of the three- 
dimensional cubic lattice. This derivation proceeds by 
diagonalizing the Hamiltonian using a Bogoliubov ap- 
proximation and subsequent expansion of the free energy 
in the superfluid velocities. 

The band dispersions and interaction parameters 
needed to evaluate the superfluid-drag coefficient are ob- 
tained numerically by means of a truncated plane-wave 
expansion of the optical lattice potentials. Using this 
numerically exact approach, instead of the tight-binding 
approximation, we obtain results that are reliable even 
in the limit of shallow optical lattices where superfluidity 
is assured for both commensurate and incommensurate 
filling. As a result of the interplay between the interac- 
tion and kinetic energies, we find a non-monotonic de- 
pendence of the superfluid-drag coefficient on the optical 
lattice depth. In contrast to Ref. [ID, where the mass ra- 
tio which maximizes the superfluid drag was found to be 
around unity for an arbitrary lattice depth, we find that 
the optimal mass ratio depends on the lattice depth. 

In addition to the two-dimensional square lattice, we 
investigate the drag in the particularly interesting three- 
beam lattices (3BL). These are two-dimensional optical 
lattices with non-separable potentials [22| . created by 
three in-plane laser beams [23[ . One special case of 3BLs 
is the triangular optical lattice, which has lately received 
attention in connection with the experimental observa- 
tion of the SF-MI transition [24j ]. 

The paper is organized as follows. In Sec. [IT] we present 
the theoretical description of two-species BECs in opti- 
cal lattices. At the same time we introduce the notations 
and conventions to be used throughout the paper. We 
discuss the two-species Bose- Hubbard model in Sec. Ill Al 
Sec. Ill Bl is devoted to the numerical derivation of the pa- 
rameters of the model, while in Sec. Ill CI we specify the 
optical-lattice potentials under investigation. 
In Sec. IIIII we generalize the derivation of an expression 
for the superfluid-drag coefficient at zero temperature, 
valid for an arbitrary lattice geometry. The results for the 
superfluid drag in two different 3BLs and the four-beam 
square lattice are presented and discussed in Sec. [IVJ Fi- 
nally, we conclude in Sec. |Vj 
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II. SYSTEM AND MODEL 

In this section, we describe the model of a two-species 
BEC in a quasi two-dimensional optical lattice. We in- 
troduce the Hamiltonian of the system (Sec. Ill Al) . out- 
line the numerical evaluation of its parameters ( Sec. Ill Bj) 
and discuss the two-dimensional optical lattices under 
consideration fSec. Ill Cj) . For convenience, we set H = 1 
throughout the paper. 



A. Two-species Bose-Hubbard model 

The single-band, two-species Bose-Hubbard model 
with short-range on-site interactions reads Q 

H = E 4aU« + \ E u a pb\ a b%b lP h a . (i) 

The operator b\ a (6 JQ ) creates (annihilates) a boson of 
component a = A, B occupying a Wannier orbital cen- 
tered at lattice site Ri (i — 1,...,N). In terms of the 
(three-dimensional) Wannier functions W a (r), 



J d 3 rW*(r-Ri) 



2m r 



V(r) 



U a p= la p / d\ \W a (r)\ 2 \Wp{r)\ 2 , 



Wair-Rj), 
(2) 



where the diagonal and off-diagonal elements of efj re- 
spectively correspond to the on-site energies and the hop- 
ping amplitudes. The potential V(r) includes the lat- 
tice potential as well as a confinement in the z-direction; 
7 Q( 3 = 2ir(m a + mp)a a p / (m a mp) is determined by the 
particle masses and the s-wave scattering length a a p. 

We consider a one-dimensional optical lattice in the 
z-direction, which provides the aforementioned confine- 
ment. Approximating the confining potential around one 
of the minima leads to a harmonic potential with mass- 
dependent frequency uj za . As long as the correspond- 
ing oscillator length l z = (raj 2(I ) -1 / 2 is small compared 
to the scattering lengths, we can assume that the two- 
body collisions are not affected by the confinement and 
Eq. @ gives the correct interaction parameters. For a 
sufficiently deep lattice in the z-direction, the Wannier 
functions can be written as 



W a (r,z) = W a (r)x (^^) 



1/4 



(3) 



where r is a two-dimensional vector. Hereafter, all the 
bold letters denote two-dimensional vectors. 

The z-dependence of U a j3 can then easily be integrated 
out and we are left with 



U a f3 = 7a/3 



m a uj za mpuj zl 3 
ir(m a uj za + mpLLi z p) 



W a (r) Wp{r) 



To switch from a real-space (lattice) description to mo- 
mentum space, we use the Fourier-transformed boson op- 
erators bi a = A -1 / 2 J2k bk a e~ lkRi - Inserting this rela- 
tion into Eq. (fTJ) leads to 

H = y^gfeqfrfeaftfca (5) 
k , a 



aft fci,...,fc4 



The band dispersion Eku and the interaction parameters 
U a p are calculated numerically as described below. 



B. Plane- wave expansion 

Using the expansion of the Bloch functions in the 
reciprocal-lattice vectors G (V denotes the system vol- 
ume) 

^(r) = ^=^Cg !G e i ( fe+G )-, (6) 

* G 

the Bloch eigenvalue problem can be recast as 

J2 <fe + G\H?\k + G')Cl G , = e ka Cl G . (7) 

G' 

Here |fe) is shorthand for y- 1 / 2 e lk r j while H" denotes 
the single-particle Hamiltonian which includes the kinetic 
energy and the two-dimensional lattice potential (but not 
the confining potential). 

Taking into account a finite number Nq of reciprocal- 
lattice vectors leads to an eigenvalue problem of finite 
dimensionality. The latter can be solved numerically to 
obtain the band dispersion and the Bloch functions [25j j . 
In our case, Nq > 100 leads to the requisite numeri- 
cal precision in the diagonalization procedure for all the 
lattices. 

To evaluate the interaction parameters U a p [cf. Eq. 
((H)], we first compute the Wannier functions 



Wa {r) = -j=^ ka {r) 



(8) 



Because for each value of k the Bloch function is only 
defined up to a phase factor, the last definition of the 
Wannier function is not unique. Therefore, one has to 
choose a gauge ^>ka —> e^^'^ka which leads to localized 
Wannier functions [H,^. The substitution [H 



$ fca (r) exp [-!lmln* fcQ (0)]f te (r) 



(9) 



(4) 



gives rise to Bloch functions which have the same phase 
at r = for each k, leading to Wannier functions which 
are sufficiently localized that intersite interactions can be 
neglected. 
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C. The optical lattices 

The optical lattices we investigate are two different 
3BLs, one with triangular geometry (3BTL) and one with 
square geometry (3BSL), and the four-beam square lat- 
tice (4BSL). Hereafter, the four-beam and three-beam 
lattice potentials are respectively labeled by the super- 
scripts 4B and 3B: 



V iB (r) 
V 3B (r) 



2 
2 



cos(2kLx) + cos(2kLy) 
cos(bi • r) + cos(62 • r) 
+ cos([6i + b 2 ] -r) 



(10) 



Here bi = ki — ki + i are the reciprocal-lattice vectors, 
where ki denotes the wave- vector of the i-th laser beam 
with magnitude = 2tt/Xl. To produce the desired lat- 
tice geometries, the lasers need to be red detuned [Vq < 
in Eq. (fTUj) ] . The magnitude of Vq is in principle depen- 
dent on the atomic properties of the Bose components. 
In this work, however, we consider systems where Vq is 
approximately equal for both species. This is not an 
unreasonable assumption since such systems are experi- 
mentally accessible (14| . 

The 4BSL is created by four laser beams of equal in- 



tensity and polarization enclosing mutual angles of tt/2 
and has a lattice period of For the 3BSL, the three 

laser beams enclose the angles tt/2, tt/2, and tt resulting 
in a lattice spacing of \ L /V2, while for the 3BTL all 
the enclosing angles are 27r/3 and the lattice spacing is 
2Al/3. The different lattice geometries are shown in Fig. 
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FIG. 1: (Color online) Optical-lattice potentials in units of 
the recoil energy Er — k 2 /2m a with their minima set equal 
to zero, (a) 3BTL, (b) 3BSL, (c) 4BSL. For these plots V = 
-En/2. 



Straightforward evaluation of the matrix elements of 
Hi in the plane-wave basis yields 



(k + G^H^k + G^) 



' 'k + G,J 2 



2m a 



In 



(k + Gi uni \V (r)\k + G htn2 ) = —[S h j 2 ±iS ni . n2 + S h ,i 2 5 nu n 2 ±ij , 
(k + G[ uni \V 3B (r)\k + G; 2i „ 2 ) = -^-(Si u i 2 ±i6 nun2 + 5i u i 2 6 nu , l2 ± 1 + Si u i 2 ± 1 6 nun2 ±i^j , 



(11) 



r 



where Gi iU = lb\ + nb 2 (/, n are integers) are the 
reciprocal-lattice vectors. The diagonal matrix elements 
are independent of the lattice potential, while the off- 
diagonal ones only depend on the potential and are equal 
for all 3BLs. 



We then neglect all the terms of higher order than bi- 
linear in the creation/annihilation operators of fc 7^ 
states. The resulting bilinear Hamiltonian can be diag- 
onalized by following the procedure outlined in Ref. l3Cl 
In terms of the Bogoliubov quasiparticle operators /3fc CT , 
one obtains 



III. DERIVATION OF THE 
SUPERFLUID-DRAG COEFFICIENT 



The derivation treated in this section generalizes 
Ref. UH details are therefore omitted. Considering the 
case where the ground state (fc = 0) is macroscopically 
occupied, the corresponding boson creation and annihi- 
lation operators commute to a very good approximation. 
Thus we can replace them by c- numbers [29| : 



{bo a bl a ) = (bt a b 0a ) + l^N 0a = (bl a b 0a ) 



boa 



u 0a 



(12) 



where the primed sum runs over all k ^ states and the 
two-branch (cr = ±) excitation spectrum is given by 
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£ ka = -L {e£ (e£ + 2F A ) + e B k (ef + 2F B ) + ayj [e£ (e£ + 2F A ) - e* (ef + 2i^)] 2 + 16F| B e£e£ 



1/2 



(14) 



r 



Here F a = n a U aa , F AB — ^Jn A n B U AB , = £fc Q — £oa, 



and 



e fe + F a ] n a = N a /N is the particle density 
of component a. It is easy to see that this spectrum 
is gapless. Since the Bogoliubov approximation makes 
the assumption of a small condensate depletion, only the 
supcrfluid phase can be described by the obtained quan- 
tities du. 

For small superfhhd velocities v a , the free energy of a 
two-component Bose-gas can be expanded as [16| 



F = F Q 



V 
2 



P>\ + pWb 



pd (v A - v B y 



(15) 



where p s a is the superfluid density of component a, while 
Fq denotes the terms independent of the superfluid ve- 
locities. At zero temperature, the free energy is equal to 
the expectation value of the Hamiltonian with no quasi- 
particles excited 



F 



T=0 



<#0> 



£k- 



EE) 



(16) 



The superfluid density of a one-component system is 
determined by its response to an externally induced su- 
perfluid velocity [HJ . To determine the drag, the corre- 
sponding two-component generalization of this procedure 
reads \\A 



el -> e%-m a v a -V k e% + 0(v 2 a ) 



(17) 



Using this transformation, the superfluid-drag coefficient 
can be found by expanding the free energy in the super- 
fluid velocities. Considering the case of parallel super- 
fluid flows, we find 



(18) 

where u denotes the direction of the superfluid flow and 
dk u = u ■ Vfc stands for the corresponding directional 
derivative. Note that the drag is always positive and 
independent of the sign of the interspecies interaction. 
The last expression holds for an arbitrary lattice. 



p A = N A m A /V. Superfluidity at all lattice depths is as- 
sured by choosing an incommensurate filling n A = n B = 
V2 [9j. The intraspecies scattering lengths are set to 
a AA — 100 ao (ao is the Bohr radius) and a BB — 65 a®. 
We choose the superfluid flows to be codirected in the 
x-direction defined in Fig. [T]for all the lattices. For the 
laser wavelength, which determines the lattice spacing, 
we choose Xl — 1064 nm. The results obtained for the 
superfluid drag are shown in Fig. [5J Figures HJb), (d), 
(f) correspond to an interspecies interaction just below 
the value where phase separation occurs [i.e., where the 



Pd/PA* 10 
4| 




144 
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IV. RESULTS AND DISCUSSION 

In the following, we present our findings for the su- 
perfluid drag in different lattice geometries, for weak 
and strong interspecies scattering. To obtain a dimen- 
sionless quantity, we normalize the superfluid drag by 



FIG. 2: (Color online) Superfluid drag for different lattices 
and interaction strengths: (a) and (b) three-beam triangular 
lattice (3BTL), (c) and (d) three-beam square lattice (3BSL), 
(e) and (f) four-beam square lattice (4BSL). The left column 
[(a), (c), (e)] corresponds to weak interspecies interactions 
a-AB = 30 ao, the right column [(b), (d), (f)] to strong ones 
aAB = 64 ao. Note the different scales for the color scheme. 
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FIG. 3: (Color online) Superfluid drag in a 3BTL as a func- 
tion of the lattice depth Vq for fixed mass ratios corresponding 
to the mixtures 87 Rb- 85 Rb (m B /m A « 1, solid), S7 Rb- 41 K 
(m B /m A « 2.2, dashed) and 87 Rb- 23 Na (m B /m A « 3.8, dot- 
ted). Component B corresponds to 87 Rb in all three cases. 
The interspecies scattering length is set to clab = 64 ao. 



excitation spectrum in Eq. (|14|) becomes imaginary]. 

Keeping the intraspecies interactions constant, the 
drag effect is enhanced significantly with increasing inter- 
species interactions. In contrast to the results of Ref. Il8l . 
however, we find that the mass ratio which maximizes the 
drag depends on the lattice depth, varying significantly 
from unity as one goes to shallower lattices (\Vq\ < I En). 



Pd -{Nm A /V)- 1 x 10 4 




For clarity, cuts of Fig. [2fb) for experimentally relevant 
mass ratios are presented in Fig. [3] 

Although the magnitude of the superfluid drag varies 
strongly with the lattice geometry, the qualitative behav- 
ior as a function of the mass ratio and the lattice depth 
is the same: for fixed mass ratio, the drag increases upon 
raising \Vq\ from zero, reaches its maximum and then 
decreases. Quantitatively, this behavior is strongly de- 
pendent on the mass ratio. 

The lattice-depth dependence of the drag can be as- 
cribed to a competition between the interaction and ki- 
netic energies. As \Vq\ increases, the interactions be- 
come stronger and the drag, which depends quadrati- 
cally on Fab, is enhanced. At the same time, the band- 
width, incorporated in pd through the derivative of the 
single-particle dispersion, decreases with increasing lat- 
tice depth. This effect is strongly dependent on the par- 
ticle masses involved. A similar behavior can be observed 
upon increasing tub for a fixed lattice depth. For a small 
mass ratio the drag increases, because it is linear in tub, 
before it starts to decay due to the decrease in kinetic 
energy. 

Comparing the 3BSL with the 4BSL, one notices that 
in the 4BSL the superfluid drag is much larger over a 
broader range of lattice depths. This is consistent with 
the fact that the hopping amplitudes are smaller and the 
SF-MI transition sets in at much shallower lattices for 
the 3BLs [23| • A possible explanation as to why the drag 
is stronger in the 3BSL than in the 3BTL can be given if 
one envisions the drag as being mediated by component 
A particles, dressed by a cloud of component B particles 
and vice versa. In this case one could argue that the 
higher coordination number of the 3BTL leads to more 
events whereby one of the particles within the cloud splits 
off the dressed particle. 

In addition to the above findings, we stress that for a 
fixed mass ratio, the lattice depth Vq which maximizes 
the drag seems to be independent of the particle densities 
n a . On the other hand, the magnitude of pd depends on 
the particle densities as shown in Fig. @] At fixed Vq, 
this dependence itself varies with the mass ratio. Com- 
paring the density dependence at Vo which maximizes the 
drag for the respective mass ratio, we find that it hardly 
changes at all. 

Finally, we note that for two components with the 
same scattering length the drag always increases upon 
increasing the particle number (data not shown). For 
two species with different scattering lengths on the other 
hand, increasing the number of the more strongly inter- 
acting particles can actually lead to a decrease in the 
superfluid drag (see Fig. 2]). 



FIG. 4: (Color online) Superfluid drag in a 3BTL (not nor- 
malized with tia) as a function of the particle densities for 
A = 85 Rb and a lattice depth Vo = —1.2 Er which maximizes 
the drag for this mass ratio (see Fig. The other param- 
eter values are the same as in Fig. [2] (b). Component B 
corresponds to the more weakly interacting species. 



V. CONCLUSIONS 

We have studied the drag between components of a 
two-species BEC in optical lattices. To assure that our 
results are valid even for shallow optical lattices, we did 
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not make use of the tight-binding approximation but used 
a numerically exact approach instead. Consequently, 
we generalized the previously derived expression for the 
supcrfluid-drag coefficient to arbitrary lattice geometries. 
To clarify the dependence of the superfluid drag on the 
lattice geometry, we have presented results for rectangu- 
lar and non-rectangular lattices with separable and non- 
separable potentials. 

We have demonstrated a non-monotonic dependence 
of the drag on the lattice depth that results from the 
competition between two effects: the drag increases with 
the interspecies interaction strength and is reduced upon 
decreasing the kinetic energy. While the qualitative be- 



havior of the drag is the same for all the lattice geometries 
studied, its quantitative properties, such as the magni- 
tude, differ from one lattice to another. 

Our study will hopefully motivate drag experiments 
with ultracold atoms in optical lattices. 
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